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Preface 



The ongoing search of extrasolar planets is one of the most attractive fields of research in 
astrophysics and astronomy. Up to now, 360 extrasolar planets have been discovered near stars 
with similar mass as the Sun. There is also discovery related to the so-called Earth-like planets. 
With regards to these discoveries, one intriguing question is whether there is relationship between 
orbit distance of the planets and their stars. Various formulas have been suggested since 1990s, and 
they suggest that there may be reason to accept quantization of distances of those planets both in 
our solar system and also in extrasolar systems as well. This book discusses this issue (Rubcic & 
Rubcic], along with other interesting issues such as protoplanetary formation of solar system 
(Pintr, prof. Perinova, & dr. Luks], precession in solar system (Pitkanen] and other topics. 

Another line of thought explored herein is the correspondence between cosmological phenomena 
and condensed matter physics, and therefore we can think that the quantization of orbit distances 
can be caused by superfluid helium quantization. This issue is explored by F. Smarandache and V. 
Christianto. Moreover, F. Smarandache also discusses possible new era of research that is 
pertaining to superluminal physics and instantaneous physics. Ion Patrascu and D. Rabounski 
discuss superluminality from their perspectives. And M. Pereira discusses his Hypergeometrical 
Universe model. 

This book is published after our previous book: Quantization in astrophyisics, Brownian motion, 
and Supersymmetry which was released about five years ago. Perhaps the ideas presented herein 
will have impact on discussions concerning quantum cosmology, which so far it cannot be observed. 
On the contrary, quantization at large scales can be observed. We hope that this volume will add a 
new chapter in our understanding of the Universe, from the viewpoint of quantization and 
discretization at large scales. 

Special thanks go to journal editors who have granted permission to reprint papers included here, 
including Chaos, Soliton, Fractals editor, Prespacetime Journal editor, Fizika editor, Progress in 
Physics editor and Apeiron editor. 



January 7 th , 2012, January 25 th 2012 
FS, VC, PP 



www.sciprint.org 
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PLANETARY ORBITS IN SOLAR AND EXTRASOLAR SYSTEMS 

ANTUN RUBCIC and JASNA RUBCIC 

Department of Physics, Faculty of Science, University of Zagreb, Bijenicka c. 32, 
10000 Zagreb, Croatia E-mail address: rubcic@phy.hr 



The analysis of orbital parameters of planets and main planetary satellites of the 
solar system, published by the present authors, suggests that the Sun’s system 
could be a “prototype” for the distribution of orbits in extrasolar planetary sys- 
tems. Owing to the recent endeavours in detecting exoplanets, it became possible 
and suitable to check this assumption. Particularly useful in this work are the mul- 
tiple extrasolar system with at least four planets. Unfortunately, there are only 
four stars satisfying this requirement. At the present time eleven stars with three 
planets have also been observed, which may also be taken into account in reaching 
reasonable assertions. Quantization of orbits in the solar system by orbital number, 
the integer n, and quantization of the product nv n (v n is the orbital velocity) by 
the spacing number, integer k, is also found in extrasolar planets. It is expected 
that new discoveries will support the present findings. 

PACS numbers: 95.10.Ce, 95.10.Fh, 95.30.-t UDC 523.2, 531.35 

Keywords: quantization of orbits, solar planets, satellites of planets, extrasolar planetary 
systems 



1. Introduction 

In our previous articles [la, b, c, d], the square law for orbits has been deduced 
by the analysis of orbital parameters of Sun’s planets and main satellites of Jupiter, 
Saturn and Uranus. The Sun’s planets are classified in two subsystems: terrestrial 
and Jovian. Therefore, there are five subsystems in the solar system, for which the 
orbital distributions follow the square law in the form 

r n = Tin 2 . ( 1 ) 

The values of n are consecutive integer numbers in a definite range and rq is the 
radius of the orbit with n = 1, dependent on the subsystem. The existing orbits, 
as an example, for terrestrial planets are distributed from n = 3 for Mercury and 
ending with n = 8 for Ceres. Similar results are obtained for other systems, as will 
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be shown later in relevant graphs and tables. In the terrestrial system of planets, 
the Earth’s moon, by hypothesis, had its primordial orbit n = 7 between Mars 
(n = 6) and Ceres (n = 8). Thus, the Moon is considered to be a planet, which 
was captured by the Earth [le] . This hypothesis is supported by the analysis of 
masses, volumes and periods of all terrestrial planets. Why the Moon at orbit 
n = 7 migrated through the orbit of Mars to become a satellite of the Earth at 
n = 5 is not clear, as well as the problem of its chemical constitution. If the Moon 
was born at orbit n = 7, then it would be expected to contain a significant amount 
of water, like both Mars and Ceres. However, the absence of volatile elements and 
water in Moon’s materials brought by astronauts suggests another origin: a giant 
impact of a body as large as Mars with the Earth. But it is hard to accept that such 
a cataclysmic process could have resulted in the Earth’s satellite with parameters 
compatible with that of the present Moon and with those of all terrestrial planets. 
Therefore, the problem of the origin of the Moon remains open. Here we treat the 
Moon as a small planet of the terrestrial group of planets. 

Physical basis for the square law (1) is a quantization of the specific angular 
momentum of planets. Details are presented in Ref. [1]. Equation (1) in extended 
form [lc,d] is given by „ 

1 77 ^ 

Tn = ^2 ' ( 2 ) 

Other relevant relations are: 



specific angular momentum 



Jn 

m n 



2gm" t 

vq k 



(3) 



1 n 3 

period T n = 2n-^GM — , (4) 

v o 
k 

velocity v n = v 0 - , (5) 

n 

where n is the orbital number, G is the universal gravitational constant, M is the 
mass of the central body, vq is the velocity constant for all subsystems in the solar 
system (close to 24 kms^ 1 ), and k is spacing number that depends on the system 
and defines the packing of orbits. In these formulae, circular orbits are assumed 
with radii equal to the semi-major axes of actual orbits. 



For a definite value of k, Eqs. 2 to 5 are simply: 



r'J 2 ~n, 


(2 7 ) 


Jn/m. n ~ n, 


(30 


T^-n, 


(T) 


nv n = const. 


(50 



It is important to point out that in a given system riv n is constant. In the solar 
system, there are five subsystems, but each with its own value of nv n . These values 
are determined by the number k. Jovian planets and satellites of Uranus have k = 1 
and almost equal values nv n . For terrestrial planets, k = 6 and similarly for other 
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systems (see Table 1). It means that physical laws are equal in planar gravitational 
systems regardless the mass m n of orbiting bodies and the mass M of the central 



TABLE 1 . Solar and extrasolar planetary systems with at least four planets. The 
masses are expressed in terms of mass of Jupiter (Mj) or Earth (Me). T is the 
period of rotation, a the semimajor axis, n the arbital number, nv n the product of 
the orbital number and velocity and k is the spacing number. 



System 


Mass 


T 

(days) 


a 

(AU) 


n 


nv n 

(km/s) 


k 


HD 160691 


c 


0.0332 Mj 


9.638 


0.09094 


i 


102.6 






d 


0.5219 Mj 


310.55 


0.921 


3 


96.8 






b 


1.676 Mj 


643.25 


1.5 


4 


101.5 






e 


1.814 Mj 


4205.8 


5.235 


7 


94.8 






(nv Il )=99±3 


4 


55 Cnc 


(e) 


0.024 Mj 


2.81705 


0.038 


1 


(146.8) 






b 


0.824 Mj 


14.65162 


0.115 


2 


170.8 






c 


0.169 Mj 


44.3446 


0.24 


3 


176.6 






f 


0.144 Mj 


260 


0.781 


5 


163.4 






d 


3.835 Mj 


5218 


5.77 


14 


168.4 






(nv n )= 170±4 


7 


GI 581 


e 


0.006104 Mj 


3.14942 


0.03 


3 


310.9 






b 


0.0492 Mj 


5.34874 


0.041 


4 


332.3 






c 


0.01686 Mj 


12.9292 


0.07 


5 


294.5 






d 


0.02231 Mj 


66.8 


0.22 


9 


322.5 






(nu„)=315±10 


13 


Ter.pl. 


Me 


0.0056 M E 


87.96 


0.387 


3 


143.6 






V 


0.815 M e 


224.70 


0.723 


4 


140.0 






E 


1 M e 


365.26 


1 


5 


148.9 






Ma 


0.107 M e 


686.98 


1.524 


6 


144.8 






Moon? 


0.012 M e 


1089 


2.07 


7 


144.8 






Ce 


0.00016 M e 


1680 


2.77 


8 


143.5 






(nv n )= 144±2 


6 


Jov.pl. 


J 


318 M e 


4333 


5.203 


2 


26.1 






S 


95 M e 


10759 


9.54 


3 


28.9 






u 


14.5 M e 


30685 


19.18 


4 


27.2 






N 


17.2 M e 


60188 


30.06 


5 


27.2 






PI 


0.002 M e 


90700 


39.44 


6 


28.4 






(m?„)=28±l 


1 


HD 10180 


b(?) 


1.35 M e 


1.17768 


0.02225 


1 


205.542 (?) 






c 


13.10 M e 


5.75979 


0.0641 


2 


242.147 






d 


11.75 M e 


16.3579 


0.1286 


3 


256.586 






e 


25.1 M e 


49.745 


0.2699 


4 


236.108 






f 


23.9 M e 


122.76 


0.4929 


5 


218.409 






g 


21.4 M e 


601.2 


1.422 


9 


231.588 






h 


64.4 M e 


2222 


3.4 


14 


233.058 






(m?„)=236=t9 


10 
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body, provided that m n are much smaller than M . Consequently, it is expected 
that the stars with their own planets must also follow the same physical laws. A 
premonition of that statement was given in Ref. [lc], but at the time of publica- 
tion, only three extrasolar systems, each with three planets, had been detected: 
PSR B1267+12, PSR1828 11 and v Andromedae. Obviously, insufficient observa- 
tional data could not present a convincible proof. However, nowadays there are 11 
detected multiple extrasolar system with thee planets (Table 2), two systems with 
four planets, one system with five planets and one system (the most recently discov- 
ered, Ref. [5] update December 2010) with seven planets. The systems with at least 
four planets (Table 1) are the best to confirm the square law for the distribution 
of orbits. This will be discussed in the next section. 



2. Analysis of observational data 

Orbital radii of terrestrial and Jovian planets and of main satellites of Jupiter, 
Saturn and Uranus are distributed according the square law (1). The circular or- 
bits are assumed with radii equal to semi-major axes. The fundamental physical 
reason is the quantization of the specific angular momentum, which in the used 
approximation is given by 



Jn/m n = ( GMn) 1 / 2 n . 

For elliptical orbits, this relation is given by [4] 

Jn/mn = [G(M + m n )a n { 1 - e 2 )] 1/2 , 

where a n is the semi-major axis and e n the eccentricity related to the n-th orbit 
(subcript is added by the present authors). For m n <C M and small values of e n , 
the approximation of circular orbits is very good. 

Using the model defined by Eqs. (1 5), or Eqs. (2' -5') for a given k = const., 
the numbers n of all orbits in a system are easily determined by the following simple 

1/3 

calculation. The values of T n are each divided by one number from a choice of 
small integer numbers (see Eq. (4')) with the aim to obtain a constant quotient for 
all orbits. 

For example: for the star 61 Vir, the periods of the planets b, c and cl (Ref. 5 
and Table 2) are: T(days) = 4.215, 38.021 and 123.01. Consequently, T 1 / 3 = 1.615, 
3.363, 4.973. One simply obtains 1.615/1 = 1.615, 3.363/2 = 1,682 and 4.973/3 = 
1.658, so the orbital numbers are n = 1,2 and 3. The resulting approximate set of 
orbital periods is Tn 3 = 1.65n, n = 1, 2, 3. 

The conclusion is that this star has three planets in successive orbits with n = 1, 
2, 3. There maybe other planets with higher n that have not been detected yet. 

However, another set of possible numbers n is obtained taking 1.615/2 = 0.808, 
3.363/4 = 0.841 and 4.973/6 = 0.829, so the orbital numbers could be n = 2, 4 

1 /3 

and 6. The resulting set of orbits is then T n = 0.83n. Other orbits would then be 
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TABLE 2. Extrasolar planetary systems with three planets. 



System 


Mass 

(Mj) 


T 

(days) 


a 

(AU) 


n 


nv n 

(km/s) 


k 


61 Vir 


b 


0.016 


4.215 


0.05020 


T 


129.6 






c 


0.0573 


38.021 


0.2175 


2 


124.4 






d 


0.072 


123.01 


0.476 


3 


126.3 






(nv n )= 127+3 


5 


Ups And 


b 


0.69 


4.617136 


0.059 


1 


139.0 






c 


1.92 


241.33 


0.832 


4 


150.0 






d 


4.13 


1278.1 


2.51 


7 


149.6 






(nv n )= 146 + 6 


6 


HD 69830 


b 


0.033 


8.667 


0.0785 


2 


197.1 






c 


0.038 


31.56 


0.186 


3 


192.4 






d 


0.058 


197 


0.63 


6 


208.8 






(nv n }= 200 + 9 


8 


GLIESE 876 


d 


0.02 


1.93785 


0.021 


2 


233.6 






c 


0.83 


30.258 


0.132 


5 


237.3 






b 


2.64 


61.067 


0.211 


6 


225.5 






(nv n )= 232 + 6 


10 


HD 40307 


b 


0.0132 


4.345 


0.047 


3 


355.8 






c 


0.0216 


9.62 


0.081 


4 


366.4 






d 


0.0288 


20.46 


0.134 


5 


356.3 






(nv ra )=360±6 


15 


PSR 1257+1 


b 


7e-05 


25.2620 


0.19 


5 


409.1 






c 


0.013 


66.5419 


0.36 


7 


412.0 






d 


0.012 


98.2114 


0.46 


8 


407.6 






(nv ra )=409±3 


17 


HD 181433 


b 


0.238 


9.3743 


0.08 


1 


92.8 






c 


0.64 


962 


1.76 


5 


99.5 






d 


0.54 


2172 


3.0 


6 


90.2 






(mv)=94±5 


4 


HD 74156 


b 


1.88 


51.65 


0.294 


1 


61.9 






d 


0.396 


336.6 


1 


2 


65.3 






c 


8.03 


2476 


3.85 


4 


67.7 






(nv n )= 65±3 


3 


HD 37124 


b 


0.64 


154.46 


0.529 


2 


74.5 






d 


0.624 


843.6 


1.64 


3 


64.5 






c 


0.683 


2295 


3.19 


4 


60.5 






(nv n )= 68 + 7 


3 


HIP 14810 


b 


3.88 


6.67386 


0.0692 


1 


112.9 






c 


1.28 


147.73 


0.545 


3 


121.7 






d 


0.57 


962 


1.89 


5 


108.0 






(nu„)=114±7 


5 


47 Uma 


b 


2.53 


1078 


2.1 


3 


63.6 






c 


0.54 


2391 


3.6 


4 


65.5 






d 


1.64 


14002 


11.6 


7 


63.1 






(nv n )= 64±2 


3 
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possible. For n = 1, the period would be Tj = 0.57 days. In Ref. [5], there are no 
stars with planets with such a small period (only the star WASP 19 has the planet 
b with the period equal to 0.789 day [5]). 

Similar calculations may be performed with semi-major axes. For the planets 
b, c and cl of the star 61 Vir, the semi-major axes (in AU) are 0.0502, 0.2175 and 
0.476 (Table 2). Then the values of a 1 / 2 are 0.224, 0.466 and 0.690. It follows that 
0.224:1= 0.224, 0.466:2= 0.233 and 0.690:3= 0.231. As expected from the previous 
considerations, the orbital numbers are n = 1,2, and 3 and col’ 1 = 0.23n. As above, 
another possible set of numbers for known orbits is n = 2, 4, 6, then a n = 0.115n, 
which for n = 1 gives cq = 0.013 AU. Note, there are no stars with such a small 
planetary orbit (only star GJ 1214 has planet b with the semi-major axis a = 0.014 
AU [5]). The number of missing planets would be three, at n = 1, 3 and 5. Which 
set of numbers is real has to be confirmed by additional observational data. This 
ambiguity is at present unavoidable. 

Another example is the star 47 Uma with three detected planets, with assigned 
orbital numbers n = 3, 4, and 7 (see Table 2). Square root of a n of the planets b, 
c and d are 1.45, 1.90 and 3.41 (AU) 1 / 2 . Following Eq. (2'), one obtains: 1.45/3 = 
0.48, 1.90/4 = 0.48 and 3.41/7 = 0.49. Consequently, occupied orbits are n = 3, 4 
and 7 and orbits 1, 2, 5 and 6 have not been detected, or may even be nonexisting. 
Again, a definitive conclusion may be obtained only by new observational data. 

In the examples above, the systems with three planets are discussed in which 
the method of the determination of orbital numbers in a planetary system is given. 

In the following analysis, the systems with at least four planets are exam- 
ined. These are the Sun’s terrestrial and Jovian planets, and the planets of stars 
HD160691, 55 Cnc, GI 581 and HD 10180. 

Table 1 shows the parameters of the planets: the masses, periods, semi-major 
axes, calculated orbital numbers n, and the products nv n with average values and 
errors, and for each system the spacing number k (last column). 

Orbital velocity is calculated using the simple formula v n = 2ira n /T n AU/day, 
in which it is assumed that a n could be taken as the radius of circular orbits. If a n 
in this formula is in units AU and T n in days, then v n = 1.0879 10 4 a n /T n kms 1 . It 
was shown that nv n is nearly constant for a particular system according to Eq. (5), 
but depends on the value of the integer k [lc] . 

The dependence of T 1 / 3 on n is illustrated in Fig. 1. 

Figure 1 are also shows the data for the satellites of Jupiter, Saturn and Uranus, 
also given in Table 3. This is done according to our statement that in all planar 
systems, the existing bodies rotate about the central large body in orbits according 
to the same physical law, and in particul satisfy the quantization of the specific 
angular momentum. It means that Jupiter with its satellites may be considered 
as a small planetary system. That similarly holds for other systems. For example, 
planetary system of the star Cnc 55 “has some basic structural attributes found in 
our solar system” [6]. It is also pointed out that the HD 10180 planetary system 
shows the regular pattern of planets’ orbits, as is also seen in the solar system [7]. 
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Fig. 1. Third roots of periods divided by chosen small integers n give the straight 
lines for all bodies in each planetary system. The ns are the orbital numbers (see 
Eq. (5')). Systems with at least four bodies are shown. 

TABLE 3. Systems of satellites of Jupiter, Saturn and Uranus. 



System 


T 

(days) 


a 

(AU) 


n 


nv n 

(km/s) 


k 


Jupiter 


Am. 


0.498 


0.00121 


2 


52.65 






Io 


1.569 


0.00282 


3 


52.01 






Eu. 


3.551 


0.00449 


4 


54.98 






Gan. 


7.155 


0.00715 


5 


54.42 






Call. 


16.689 


0.0126 


6 


49.23 






{nv n )= 52.7A1.6 


T 


Saturn 


Jan. 


0.693 


0.00101 


6 


95.11 






Mini. 


0.942 


0.00124 


7 


100.24 






Enc. 


1.370 


0.00159 


8 


101.07 






Teth. 


1.888 


0.00197 


9 


102.16 






Dione 


2.737 


0.00252 


10 


100.27 






Rhea 


4.518 


0.00352 


11 


93.31 






(nv n ) =98. 69 ±3.0 


4 


Uranus 


Ariel 


2.520 


0.00128 


5 


27.55 






Umb. 


4.144 


0.00178 


6 


28.04 






Tit. 


8.706 


0.00291 


7 


25.49 






Ober. 


13.463 


0.0039 


8 


25.21 






Mir. 


1.414 


0.000865 


4 


26.62 






Puck 


0.672 


0.000575 


3 


24.62 






(nv n )= 26.26A1.2 


T 
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The dependence of nv n on k is shown in Fig. 2. Straight line obtained by linear 
regression is 

nv n = (23.64 ± 0.32)fc + 3.76 ± 0.87 kms" 1 (6) 

Since nv n = kv o 

v 0 = (23.64 ±0.32) kms -1 . (7) 



Note that nv n for the planet 55 Cnc, e is considerably smaller than those of the 
other planets with n = 2, 3, 5 and 14. The reason for that remains unknown, and 
in Table 1 both e and nv\ are given in parentheses and are not included in mean 
value (nv n ). 




Fig. 2. Products nv n of orbital number n with orbital velocities v n for all systems 
represented in Fig. 1 are shown in steps defined by the spacing number k. 

The number k defines the spacing of orbits in a system. It is interesting to 
point out that Jovian planets and satellites of Uranus both have k = 1. It means 
that orbital velocities decrease with n equaly in both systems. Thus, for the n — 5 
planet Neptune and Uranian satellite Ariel have the same orbital velocities. This is 
wonderful having in mind that the planar systems of Sun and Uranus are mutually 
nearly orthogonal. It is also impressive that planets of the star HD 160691 have 
k = 4 as the satellites of Saturn. Deviation of Vq in Eq. (6) is less than 2%, owing 
to the analysis of four or more planets per system. 
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However, 11 systems with only 3 planets per system have greater dissipation of 
nv n , as can be seen in Table 2, where maximum errors are included. Nevertheless, 
the mean nv n values may be distributed so that they are close to a straight line 
defined by Eq. (6), as is illustrated in Fig. 3. 




Fig. 3. Dependence of nv n on k for all systems shown in Fig. 2 and for all systems 
with three planets. 

The straight line determined by all considered systems listed in Tables 1 - 3 is 

nv n = (23.76 ± 0.16)fc + (2.48 ± 0.69) kms -1 . 

Therefore, Vq = (23.8 ± 0.2) kms -1 . In spite of the satisfactory description of 
the dependence of nv n on fc, the systems with 3 planets cannot safely confirm 
the quantization of nv n with the above step of Vq. Additional stars with more 
planets should be decisive for a final conclusion. Hopefully, advanced technique of 
observations of extrasolar planetary systems will help in reaching a definite solution. 



3. Conclusion 



We applied our model of quantization of orbits in the solar system on newly 
discovered extrasolar planetary systems. We confirmed that the square law for the 
distribution of orbits deduced from observational orbital parameters in solar system 
(Eq. (1) and/or Eq. (2)) also holds for the extrasolar planetary systems. Both 
integers, the orbital number n and the spacing number k , can easily be determined 





QUANTITATION AND DISCRETIZATION AT LARGE SCALES 



10 



from observed periods and semi-major axes of planets in systems considered. All 
this may be useful in the classification of orbits. Thus third root of the period 
divided by some small integer number n needs to be nearly constant for all planets 
in the system. Then, n is the number of orbit. We emphasize that the square 
law of orbits defines only the architecture of the planetary system, but details 
can only be determined by using observational parameters of some real objects 
belonging to the system considered. When the possible set of orbits is defined on 
the basis of occupied orbits, then one can anticipate which empty orbits could 
contain unobserved planets. 

For example, the 55 Cnc planets e, b, c, f, and d are, according to our analysis, 
located at orbits 1, 2, 3, 5 and 14. Eight orbits at n = 4, 6, ..13 are empty. The first 
thought is that at orbit 4 could be a small yet undetected planet. Moreover, the 
authors in Ref. [6] presume that in the gap between periods of 260 days to 13 yr 
several planets could exist and probably maintain dynamical stability. 

In the HD 10180 planetary system, the first five orbits are occupied by the 
planets b, c, d, e, and f. The last planet h with relatively large mass is at the orbit 
14. But at orbit 9 there is the planet g. Similarity with 55 Cnc system is impressive. 

The procedure outlined above has also been applied to the origin of the Moon. 
Namely, in the terrestrial planets, the orbit 7 is empty and is located between Mars 
and Ceres. We have put forward a hypothesis that the Moon originated at that 
orbit and later on migrated to be captured by the Earth [le] . The argument for 
such an assertion is that definite mass and volume of the Moon are expected when 
compared with the same quantities of all terrestrial planets. 

We hope that determination of possible orbits according to square law could be 
a guide in the search for extrasolar planetary systems. 
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STAZE PLANETA U SUNCEVOM I IZVANSUNCEVIM SUSTAVIMA 

Analiza parametra putanja planeta i glavnih planetarnih satelita u suncevom sus- 
tavu, objavljena u nasim prijasnjim radovima, upucuje na to da bi suncev sustav 
mogao biti prototip i za planetarne sustave zvijezda slicnih Suncu. Zahvaljujuci 
novijim rezultatima u detekciji izvan suncevih planeta (exoplaneta) omogucena je 
provjera ove pretpostavke.U tu svrhu su najpogodniji sustavi sa cetiri i vise planeta, 
ali nazalost takvih sustava je otkriveno samo nekoliko. Veci broj sustava ima samo 
tri planete, ali i njhova analiza daje potporu gornjoj pretpostavci iako sa manjom 
vjerodostojnoscu. Kvantizacija putanja u suncevom sustavu s cijelim brojem n, i 
kvantizacija prostornosti (pakiranja) putanja sa cijelim brojem k , vodi na relaciju 
nv n = kv o, gdje je v n brzina planete na putanji, a vq je konstantna brzina za sve 
sustave. Ove velicine mogu se odrediti i u izvansuncevim sustavima. Ocekujemo da 
ce nova otkrica exoplaneta potvrditi nasa dosadasnja saznanja. 
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APPENDIX TO“PLANETARY ORBITS IN SOLAR AND EXTRASOLAR SYSTEMS“ 



1. Orbital velocity of electron in H atom (and ionized atom He + and Li 2+ ) 

In old quantum mechanical approach to hydrogen-like atoms Bohr used the third Newton's 
law 

m e vj !r ne = (1/ 47i£ 0 )Ze 2 / rj (1) 

and quantized angular momentum, which Bohr found by his great intuition and trial and error 
method. 

m e V ne r ne = nk / 2X (2) 

In these equations m e is electron mass, v ne orbital speed of electron at n-th orbit whose radius 
is r ne , £ 0 is permitivity of vacuum, e electronic elementary charge, Z is number of elementary 
charge in atomic nucleus, and h is the Planck's constant. Equations (1) and (2) define orbital 
radius Eq.(3) and orbital speed Eq.(4) 

hn 2 



27tZmca 



(3) 



and 



nv, u = Zac 



(4) 



In Eq.(3) and (4) a=(l/47t£ 0 ) 27teVhc is the fine structure constant. 

2. In gravitational systems a similar approach to that presented by equations (1) to (4) is 
following: Third Newton's law for rotating body of mass m n at n-th orbit, of radius r n with 

in v ^ Mm 

speed v n around central body of mass M is " ” = G — y~ (5) 

r r 

n n 

G is universal constant of gravitation. 

A second necessary equation in order to solve unknown r n and v n is based on 
observational data and the square law for orbital radius of orbits in the solar system 
has been established in the form r n =rirr (6) 

As a consequence angular momentum is given by linear function on n, analogous to Eq.(2) 



t j. ,, , „M 1 v. 

It follows that the speed v n is v n = G — = G = — 

V r n V r i n n 



or 



r M 

nv n = G — =vi 



(7) 



The speed v n depends on the system considered and for several known systems it was found 
that vi = kv Q where v 0 =close to 24 kms l . Therefore, one may write 



nv„=kv 0 



( 8 ) 



Eq.(4) shows that by increasing Z the set of orbits are more densely packed and speeds v n 
increase too. 

Similarly, Eq(8) gives a growth of speeds as integer number k increases. 

It was also shown that v n =ac/27i;f 



or 



ac 

V " = T 



(9) 
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where 1/F=(0.01073±0.0013)k +(0.00188+0.00043). For k=l Vl =v o =(0.01 073+0.00 18)ac 
=27.6 kms' 1 which is valid for Jovian planets i.e. v n =27.6/n. For terrestrial planets k=6 and 
Vi=0.01073k+0.00188=145.1 kms" 1 . If the intercept is neglected in fitting procedure the resulting 
equation is 1/F=0.0 1118 k± 0.000169. For Jovian planets follows V!=24.48 kms" 1 and for terrestrial 
planets vi= 146.9 kms" 1 . However, using the tabulated values for periods and semi-major axes 
with at least four planets for multiple planetary systems the valus nv n are described by 
nv n =(23.64+0.32)k +(3.76+0.87) kms" 1 . 

where k is the integral number in a range from 1 for the Jovian system of planets up to 17 for the 
system of pulsar PSR 1257+12. Therefore, the constant 144 kms 1 is not a fundamental constant 
concerning the planetary systems of the stars similar to solar system. 

For example, planets of 55 Cnc have an average value of nv n close to 170kms"' with k=7. 

Warning: in November 2010. planet denoted by Cnc-e had period P(e) 2.81795 days and semi-major 
axis a(e) 0.038 AU., )(data are taken from Jean Schneider: Extrasolar Planets Encyclopaedia, shortly 
EPE). These data were used in submitted paper in 2010., and nv n was 146 kms" 1 . But, in December 
2011. in (EPE) one finds P)e)= 0.73054 days and semimajor axis a(e)=0.0156 AU giving nv n 230 
kms 1 . We concluded that data cited are dubious Our analysis suggests that P(e) and a(e) 
should be close to 1.807 days and 0.0285 (AU), respectively. Then, for planet Cnc e 
nv n =171.5 kms" 1 in agreement with other four planets (see Table 1) Numbers n for planets are 
1,2, 3, 5 and 14. 

Recently, the fourth planet of Gliese 876 is detected and data for all planets are listed in 
EPE. The set of orbits are numbered with 2 ,5, 6 and 8. Then the values of nv n in kms" 1 are : 
2v2=233.6, 5v5=234.3, 6v6=222.5, amd 8vs=234.1. The mean value for nv n =(231. 1+4.1) 
kms" 1 and k =10. However, this is not included in Table 1. 

In conclusion: Eq.(4) shows that fundamental speed of electron is occ=3 10 5 /137=2189 
kms" 1 . Number Z determines the packing of permissible orbits, which for all systems follows the 
square law for orbital radii. 

Eq.(8), analogous to Eq.(4), assumes the fundamental constant v G common for all 
gravitational planar systems, while integer k determines the spacing of orbits. The constant 
v o =0Cc/F=24 kms" 1 . Factor F is deduced from observational astronomical parameters and 
presently has not theoretical explanation. Consequently the speed 144 kms 1 =kv Q with k=6 is 
one of possible speed among others with various values of k in an interval from k =1 to 17, 
which is found up to now. 

The speed v G is proportional to ac and if ac is multiplied by (3/2)a then (3/2)cor =23.98 
kms" 1 is close to v Q . Obviously, this little game with numbers is not to be taken seriously, but 
could perhaps be interesting to some teoreticians in studies a connection of electromagnetism 
with gravitation. 
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Areal velocities of planets and their comparison 

Pavel Pintr, Vlasta Perinova*, Antonm Luks 
•Joint Laboratory of Optics, Palacky University, RCPTM, 

17. listopadu 12, 771 46 Olomouc, Czech Republic 



Abstract 

We have utilized the fact that the areal velocity of a planet is directly 
proportional to the appropriate number of the planet, while its distance is 
directly proportional to the square of this number. We have confirmed a 
previous proposal of the quantization of the planetary orbits, but with the 
first possible orbit of a planet in the solar system identical only to an order 
of magnitude. Using this method, we have treated moons of two planets 
and one extrasolar system. We have investigated a successive numbering 
and suggested a Schmidt-like formula in the planets and the Jovian moons. 

PACS number: 96.10.-pi 

Keywords: solar system, distances, planets and satellites 



1 Introduction 

The formation of the solar system and its development are well described in [1], 
The existing theories presume the age of the solar system as 4.5 billion years 
and that the entire system was created approximately 100 million years after the 
formation of the Sun. Despite of this, some of the chronological events of the 
formation of the system still remain unknown to us. 

Comparing with the young T Tauri stars, we can say that the Sun formed 
in the center of a protoplanetary disk with the dimensions of approximately 
1000 AU. The planets formed in the first 10 million years after the formation of 
the protoplanetary disk. The development of the solar system was terminated 
approximately 90 million years after the formation of the protoplanetary disk. 

Very interesting papers have been devoted to the mechanism in protoplanetary 
disks [2, 3]. Turbulent processes have been described in nascent protoplanetary 
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nebulas. These topics have been the subject of many papers. The beginning of 
modern theories dates since Kuiper [4], who has shown that the protoplanetary 
nebulas would have to be more massive than the algebraic sum of masses of all 
planets. 

We can divide papers describing the distribution of distances into several 
categories. Many empirical formulas describe the distances on the condition of 
suitable numbering of planets [5, 6]. Dubrulle and Graner [7, 8] have shown that 
using the rotational symmetry and the scale invariance, we can derive a geometric 
progression for any model system in the form 

r n = r 0 K n , (1) 

where n is an integral number, r'o is an initial distance, and K is a constant that 
determines the distribution of distances in the system. A successive numbering 
of planets is assumed as is respectable to such an impressive formula. Krot has 
created an evolutionary model of the rotating and gravitating spherical body 
[9]. He has remembered that with the aid of specific angular momentum of 
protoplanets, Schmidt derived the square root of radius R n of the orbit for the 
nth protoplanet [10], 

\Jlln = a + bn, (2) 

where a and b are constants. Then he has generalized the Schmidt law for the 
solar system leaving (2), a mere linear approximation. 

The first quantum formulas are comparable in complexity with (2). Agnese 
and Festa [11, 12] have described the distances of planets in the solar system 
as a gravitational atom [13] using the famous Bohr-Sommerfeld rules. The suc- 
cessive numbering is possible only for the terrestrial planets. The other planets 
are numbered “suitably”, i.e., so that a best fit is achieved. They have shown 
that this description can be applied to extrasolar systems. They have proposed 
a gravitational constant in conformity with the clue to the unification of gravita- 
tion and particle physics [14], The hypothesis of a fundamental orbital distance 
0.055 AU has been a very interesting result [11, 12, 15, 16]. A derivation of 
the Schrodinger equation [17, 18] from the Newton mechanics has inspired many 
variations of quantum description [13, 19, 20, 21], We can find solutions of the 
Schrodinger equation in [21, 22], which lead to possible discrete orbits by means 
of the quantum averaging. The distances of planets obtained in such a way ex- 
hibit a dependence on the main and orbital quantum numbers. The probability 
densities have been derived for each orbit and the number of possible orbits in 
the solar system has been reduced [22], 

Till now 360 extrasolar planets have been discovered near stars with similar 
mass as the Sun. Every day we observe new extrasolar planets or protoplanetary 
disks. Theories of migrating planets suppose that, if two high mass planets form 
near each other, both the planets will change orbits around the star and also the 
collisions with next big bodies will change orbits of planets in a young planetary 
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system [23, 24, 25]. According to these theories, the predictions of orbits are very 
problematic. 

In section 2, we will expound the method under application and use it for the 
planets in the outer part of the solar system. In section 3, we will apply it to the 
systems of moons around Jupiter and Uranus. In section 4, we will consider the 
extrasolar system HD 10 180 [26]. 



2 Correlation of areal velocities 



Agnese and Festa [11, 12] have invented allowed planetary orbits with the major 
semi-axes and excentricities 



CL n Q.‘\Yl , E n [ 




( 3 ) 



respectively, where a\ means a possible first orbit of a planet, n is a principal 
number and l is an azimuthal number, l — 1, . . . , n, n — 1, 2, 3, . . . , oo. Assuming 
/ = n (circular orbits), they describe the distribution of planetary distances in 
words that we formalize as 

rip) = hn( P ), (4) 

where p = Mercury, Venus, Earth, Mars, (Ceres,) Jupiter, Saturn, Uranus, Nep- 
tune, Pluto and n( Mercury) = 3, n(Venus) = 4, n(Earth) = 5, n(Mars) = 6, 
(n(Ceres) = 8,) n(Jupiter) = 11, n( Sat urn) = 15, n( Ur anus) = 21, n(Neptune) = 
26, n(Pluto) = 30. It can be seen that the inner planets are successively num- 
bered and the outer planets are rather numbered with a step of 5. For the first 
orbit it holds that 

GM< Sun ) 

°i = — T2~’ 5 

OigC Z 

where G is the gravitational constant, c is the speed of light, a y is a gravitational 
structure constant, 1 /a g = 2113 ± 15 [11], and M (Sun) is the mass of the Sun 
in application to the first possible orbit in the solar system. We remark that 
the gravitational structure constant, whose value was calculated from data of the 
solar system, has been tested against extrasolar planets and provided an orbit 
hi = 0.055 AU [11, 12, 15, 16]. This description has shown a very interesting 
connection between a model of the solar system and the hydrogen atom. In 
this paper, we will show that it is also possible to use quantum physics for the 
determination of the distribution of orbits in planetary systems. 

Let us consider the solar system, where we implement a simplification that 
the orbits of planets are circular and the positions of planetary orbits are in one 
plane. We can define a circular planetary model of the solar system, which looks 
like the model of hydrogen atom from the “old quantum theory”. In the old 
quantum theory, it was only possible to explain the structure of the hydrogen 
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atom or an ionized atom with a single electron. The absorption or emission lines 
for spectroscopy were obtained in terms of the energy differences of the electron 
on various orbits. For these orbits it holds that (cf. Bohr’s model of hydrogen 
atom) 

m e v n r n = nfi, (6) 

where m e is the mass of electron, v n is the velocity of electron revolving around the 
nucleus, f n is the distance of electron to the nucleus, and ft is the reduced Planck 
constant. These orbits may not be occupied simply by more and more electrons, 
because the Coulombic interaction between these particles is not negligible. This 
quantization of orbits can be formally generalized to macroscopic bodies and the 
velocities related to the gravity. Agnese and Festa [11] have modified the relation 

(6) to the form 

m p v n f n = n(h + m p cA), (7) 

where m p is the mass of planet, v n is the orbital velocity of planet, f n is its 
distance to the central star, c is the speed of light and A is a fundamental length. 

With respect to the weak equivalence principle in the case of circular orbits, 
the velocities will only depend on the gravitational potential in the distance r to 
the Sun (the central body). On neglecting ft on the right-hand side of the relation 

(7) , the independence of allowed orbits of the mass m p is obtained. In contrast 
to the electrons in the atom, the orbits around the Sun can be occupied by more 
than one macroscopic body as far as the gravitational interactions between them 
can be neglected. This ad hoc hypothesis explains the regularity of the planetary 
orbits on the given, maybe too generous assumption. 

Let us study the simplified model of the solar system and let us address to 
the following consideration, which comes out of Kepler’s second law: Areas which 
are swept out by the radius vector of planet in equal time intervals are equal, so 
the elementary area swept out by the radius vector of planet in the aphelium in 
the time df is the same as the elementary area swept out by the radius vector of 
planet in the perihelium in the time df. For the area which is swept out by the 
radius vector of planet in the circular model, Kepler’s second law is valid as well. 
For the areal velocities of planets, w(p), it holds that 

2 w(p) = v(p)r(p). (8) 

Let us compare the areal velocities of planets for the outer part of the solar system, 
with the allowed areal velocities w n , 2 w n = v n f n , which will be appropriately 
defined. 

From Table 1, we can substitute a formula v(p)r(p) = n(p)K(p ) by a new 
formula v n f n = nA"7 ipprox ), where n = n(p) and _ft'( approx ) is a constant, viz., an 
approximate value of K (p) , or 



v(p)r(p) « n(p) jf( approx ), 



( 9 ) 
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p 


v(p)r(p) 


n(p ) 


K(p) 


A'(approx) 


n(p)iC (approx) 


Jupiter 


1.02 x 10 16 


10 


1.02 x 10 15 


1.00 x 10 15 


1.00 x 10 16 


Saturn 


1.38 x 10 16 


14 


9.87 x 10 15 


1.00 x 10 15 


1.40 x 10 16 


Uranus 


1.95 x 10 16 


20 


9.75 x 10 15 


1.00 x 10 15 


2.00 x 10 16 


Neptune 


2.45 x 10 16 


25 


9.80 x 10 15 


1.00 x 10 15 


2.50 x 10 16 



Table 1: Parameters K (p) , A'( a pp tox ) anc [ n(jp) for the outer part of the solar 
system. 
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Figure 1: Comparison of real data v(jp)r(jp) (x) with the formula n(p)K ' ( approx ) 
(o) for the outer parts of the solar system. 
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where p = Jupiter, Saturn, Uranus, Neptune. A comparison of real data v(p)r(p) 
with the approximate formula n(p)Ad approx ) is illustrated in Figure 1. 

We find v n , f n such that they fulhl the relation 

v n r n = niF( approx ) (10) 



and Newton’s gravitational law 



lGM( Sun ) 



v n = 



( 11 ) 



We arrive at 

f n = apn 2 , (12) 

where 

\J{ (approx) 12 

ai = ^ r , q = 0.052 AU. (13) 

GM( S un) v ' 

We will show that the planets can be numbered successively using n(J) to 
n(J ) + 3, where J stands for the planet Jupiter, unlike 10, 14, 20, 25 in Table 1. 
Traditionally, we use the least squares method. We can determine the number 
n(J) and a constant K such that 

[1.02 - n(J)K] 2 + [1.38 - (n(J) + 1 )K} 2 + [1.95 - (n(J) + 2 )K] 2 

+ [2.45 — (n(J) + 3 )K] 2 = min. (14) 

This happens for n(J) = 2, K = 0.4857 x 10 16 . We have arrived at a Schmidt-like 
formula. 

The formula (12) is in accordance with the papers of Agnese and Festa [11, 12], 
but we do not use the gravitational structure constant for the definition of a 
possible first orbit. This is the main point in our considerations. 

Now we can introduce a new parameter p\ of the system in the formula (13), 
which we call the length density of orbits, 



G 

j ( approx) ‘ 2 



(15) 



in units kgm 1 . This new parameter can be used for the classification of extrasolar 
systems. 



3 Systems of moons around planets 

The procedure which we derived above, is valid also for systems of moons around 
planets. For a system of moons, the formula (12) is valid, where for the first orbit 
it holds that 

1 

p z M( planet ) ’ 



a i = 



(16) 
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p 


v(p)r(p ) 


n(p ) 


K{p) 


A' (approx) 


n (p) A'( a PP rox ) 


Io 


7.31 x 10 12 


7 


1.04 x 10 12 


1.00 x 10 12 


7.00 x 10 12 


Europa 


9.21 x 10 12 


9 


1.02 x 10 12 


1.00 x 10 12 


9.00 x 10 12 


Ganymedes 


1.16 x 10 13 


12 


9.67 x 10 11 


1.00 x 10 12 


1.20 x 10 13 


Callisto 


1.50 x 10 13 


15 


1.00 x 10 12 


1.00 x 10 12 


1.50 x 10 13 



Table 2: Parameters K (p ) , A'( a pp rox ) and n(jp) for the Jovian system of moons. 
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Figure 2: Comparison of real data v(p)r(p ) (x) with the formula n(p)K' P pprox ) 
(o) for the Jovian system of moons. 



where M (planet ) is a mass of a planet that moons revolve around. 

For the Jovian system, it holds that M( Jupiter ) = 1.9 x 10 2 ' kg, ai = 7890.79 
km, /i( a P prox ) = 1.00 x 10 12 m 2 s 1 . In Table 2, the parameters for the Jovian 
system of moons can be found. A comparison of real data v(p)r(p) with the 
approximate formula n(p)A'( approx ) is illustrated in Figure 2. 

We will show that the moons can be labelled with successive numbers n(I) to 
n(I) + 3, where / stands for the moon Io, unlike 7, 9, 12, 15 in Table 2. Again, 
we use the least squares method. We can find the number n(I) and a constant 
K such that 

[7.31 - n(I)K} 2 + [9.21 - (n(I) + 1)A'] 2 + [11.6 - (n(I) + 2 )K] 2 

+ [15.0 — (n(I) + 3)A'] 2 = min. (17) 

This takes place for n(I) = 3, K = 2.40 x 10 12 . We have indicated a Schmidt-like 
formula. 
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Figure 3: Comparison of real data v(jp)r(jp) (x) with the formula n(p) Ah a,pprox ) 
(o) for the Uranian system of moons. 



For the Uranian system it is valid that jvU Uranus ) = 8.7 x 10 25 kg, cp = 1723.28 
km, A'UPProx) = i go x 10 11 m 2 s 1 . In Table 3, the parameters for the Uranian 
system of moons can be found. A comparison of real data v(p)r(p) with the 
approximate formula n(p)Ah approx ) is illustrated in Figure 3. 



p 


v(p)r(p) 


n(p) 


Kip ) 


j^(approx) 


n(p) Ah approx ) 


Miranda 


8.68 x 10 11 


9 


9.64 x 10 11 


1.00 x 10 11 


9.00 x 10 11 


Ariel 


1.05 x 10 12 


11 


9.57 x 10 11 


1.00 x 10 11 


1.10 x 10 12 


Umbriel 


1.24 x 10 12 


12 


1.04 x 10 11 


1.00 x 10 11 


1.20 x 10 12 


Titania 


1.59 x 10 12 


16 


9.94 x 10 11 


1.00 x 10 11 


1.60 x 10 12 


Oberon 


1.84 x 10 12 


18 


1.02 x 10 11 


1.00 x 10 11 


1.80 x 10 12 



Table 3: Parameters K(p), A'( approx ) and n(jp) for the Uranian system of moons. 



4 Extrasolar system HD 10180 

We can also apply our consideration to extrasolar systems. For such systems it 
holds that 

[A"( a PP r0X )] 2 
Ql = GMW ’ 

where jbC Star ) is a mass of a central star. Here M^ Star )=1.06 ± 0.051bC Sur A We 
calculate further orbits according to the formula (12). 
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The planetary system HD10180 was introduced in [26]. It is the most explored 
extrasolar system with 7 planets. This extrasolar system was examined with the 
aid of measurements of the radial velocities of the system HARPS and that is 
why we selected it for our considerations. 

The planet b is at the distance 0.02226 AU to the central star with the mass 
1.4 times more than the mass of the Earth, the planet c is at the distance 0.0641 
AU to the central star with the mass 13.16 times more than the mass of the 
Earth, the planet d is at the distance 0.1286 AU to the central star with the 
mass 11.91 of the mass of the Earth, the planet e is at the distance 0.2695 AU 
to the central star with the mass 25.3 of the mass of the Earth, the planet f is 
at the distance 0.4923 AU to the central star with the mass 23.5 of the mass of 
the Earth, the planet g is at the distance 1.422 AU to the central star with the 
mass 21.3 of the mass of the Earth, the planet h is at the distance 3.4 AU to the 
central star with the mass 65.2 of the mass of the Earth. 

The formulas (3) and (12), with hi = 0.055 AU do not provide an appropriate 
allowed orbit, because the distance of the planet b to the central star r(b), is 
much nearer than the radius of a possible first orbit a For the system HD10180, 
it holds that iP( approx ) = 1.00 x 10 14 m 2 s _1 , with the first orbit hi = 0.000484 
AU and the length density of orbits pi = 6.67 x 10~ 39 kgm -1 . If we compare 
the length density of orbits with the solar system, the system HD10180 has 100 
times denser orbits than the solar system. Therefore, the architecture HD10180 
is much nearer than the solar system. The extrasolar system HD10180 meets the 
formula (12), if we apply the correct first distance hi. In Table 4, the parameters 
for the system HD10180 are arranged. A comparison of real data v(p)r(p) with 
the approximate formula n(p)Ad appr °U i s illustrated in Figure 4. 



p 


v(p)r(p) 


n(p) 


K{p) 


(approx) 


n (p)K ( approx ) 


b 


6.86 


X 


10 14 


7 


9.80 


X 


10 13 


1.00 


X 


10 14 


7.00 


X 


10 14 


c 


1.15 


X 


10 15 


12 


9.59 


X 


10 13 


1.00 


X 


10 14 


1.20 


X 


10 15 


d 


1.63 


X 


10 15 


16 


1.02 


X 


10 14 


1.00 


X 


10 14 


1.60 


X 


10 15 


e 


2.36 


X 


10 15 


24 


9.84 


X 


10 13 


1.00 


X 


10 14 


2.40 


X 


10 15 


f 


3.19 


X 


10 15 


32 


9.98 


X 


10 13 


1.00 


X 


10 14 


3.20 


X 


10 15 


g 


5.42 


X 


10 15 


54 


1.00 


X 


10 14 


1.00 


X 


10 14 


5.40 


X 


10 15 


h 


8.38 


X 


10 15 


84 


9.98 


X 


10 13 


1.00 


X 


10 14 


8.40 


X 


10 15 



Table 4: Parameters K(p), if( approx ) and n{p) for the system HD10180. 
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Figure 4: Comparison of real data v(jp)r(jp) (x) with the formula n(p)K ' ( approx ) 
(o) for the system HD 10 180. 

5 Conclusion 

Basing on the numerical agreement between calculation and real data, we have 
found suitable numbers of the planets and a proportionality constant of their 
areal velocities to this numbers. We have derived that the distance of the planet 
to the central star is directly proportional to the square of this number. In this 
way we have obtained the first possible orbit of a planet at the distance 0.052 
AU in the solar system and 0.000484 AU in the system HD10180. Analogously, 
we have got the first possible orbit of a moon at the distance 7890.79 km in the 
Jovian system of moons and 1723.28 km in the Uranian moon system. 

We conclude that the distances of the planets and moons in gravitational 
systems can be obtained as follows: 

Areal velocities of planets relate to integral numbers, viz., suitable numbers of 
the planets. These velocities are directly proportional to the appropriate numbers 
of the planets with a proportionality constant A^ approx ) . 

Distances of the planets in the gravitational system are directly proportional 
to the squares of the numbers of the planets and the proportionality constant, 
the radius of a possible first orbit d\ depends on the parameter A'( a pp rox ). 
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Distribution of distances in the solar system 
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Abstract 

The recently published application of a diffusion equation to prediction of distances 
of planets in the solar system has been identified as a two-dimensional Coulomb prob- 
lem. A different assignment of quantum numbers in the solar system has been proposed. 
This method has been applied to the moons of Jupiter on rescaling. 

PACS number: 96.35 
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1 Introduction 

The 20th century is held for the golden age of the astronomy and astrophysics, when many 
persistent questions were solved and the human view of the universe changed radically. In 
spite of this, at the beginning of the 21st century, one cannot find satisfactory answers to 
some questions our ancestors posed as early as in the 16th century. For instance, Kepler 
looked for a universal law, in his Mysterium cosmographicum, to explain the planetary 
distances in the solar system. Nowadays, when discoveries of other planetary systems occur, 
such a law could explain the distances of their planets. 

In 1766 Titius formulated the law, which described distances of the bodies in the solar 
system, and it even predicted new bodies at certain distances from the Sun [1], Actually its 
being criticized led to the discovery of the remaining planets and new bodies - asteroids - 
in the solar system. It was the first, controversial, description of the distances of the bodies 
in this planetary system. But hardly any physical explanation has thus far been given. Is 
it a mere extravagance, or does this law have some deep physical content? May the planets 
around stars originate at definite distances? 

Quest of the answer developed into invention of new empirical formulae, which describe, 
with higher or lower accuracy, the distances of the bodies in the solar system. For instance 
Armelini’s empirical formula has the form 

TnA = 1.53", (1) 

where n assumes the values: Mercury —2, Venus —1, Earth 0, Mars 1, asteroid Vesta 2, 
asteroid Camilla 3, Jupiter 4, Saturn 5, asteroid Chiron 6, Uranus 7, Neptune 8, and Pluto 
9. 

°*Corresponding author. 

E-mail address: perinova@prfnw.upol.cz (V. Perinova) 




QUANTIZATION AND DISCRETIZATION AT LARGE SCALES 



In 1938 Mohorovicic invented an empirical formula [2], which describes the distances of 
planets and comets with high accuracy, and it also predicts an asteroid belt between Mars 
and Jupiter. Mohorovicic’s law says that the distances of the inner parts of the solar system 
increase in a sublinear manner and those of the outer parts of this system increase in a 
superlinear manner. In the paper [3] we have modified this law such that it satisfies also 
other planetary systems and those of the moons of the giant planets. 

Interesting is the empirical formula, which is similar to the laws of quantum mechanics 



[ 4 ] 



1 

2 



(m 2 + n 2 )r 0 , 



( 2 ) 



where m are natural numbers, n — 0, 1, . . . , m and r 0 = 0.387 AU. The Bohr-Sommerfeld rule 
of (allowed) orbits for electrons in the electric fields of the nuclei of various atoms resemble 
the distribution of planetary distances, but do not let us forget that this rule describes 
bodies (electrons), which all have the same inertial mass and the same electric charge, which 
replaces a gravitational mass here. To obtain a distribution of the planetary distances, one 
either replaces different planetary masses by their mean mass, or makes the quantum of 
action depend on the actual mass. 

Agnese and Festa described the solar system like a gravitational atom [5]. They utilized 
a quantum law for the hydrogen atom, which they applied to description of major semi-axes 
of allowed (discretized) elliptical orbits of the bodies of the planetary system 



r n af = run 2 , 



( 3 ) 



where n are natural numbers and rq is the Bohr radius of the planetary system, which is 



T\ = 



GM 



ale 2 ’ 



( 4 ) 



where G is the gravitational constant, M the mass of the central body, c the vacuum speed 
of light and a g is a gravitational structure constant, which has the property 4- = 2113 ± 15. 
Agnese and Festa have shown that this description of distances satisfies also the planetary 
system v Andromedae [6] and other stellar systems alike on substituting the mass of the 
appropriate central star for the mass M. A study which elaborates on such ideas has been 
presented in [7]. 

Recently, the significance of the Titius-Bode law has been evaluated both by generating 
random planetary systems [8] and by the help of methods of the modern statistical analysis 
[9]. In the papers [10, 11] the authors point out quantum features also on large scales, namely 
discrete values of distances of possible planets and galaxies. 

In quantum mechanics one utilizes Schrodinger’s equation for the description of a physical 
system. In the paper [12], the stochastic mechanics is constructed, i. e., the Schrodinger 
equation is obtained as a classical diffusion equation by the help of the hypothesis that any 
particle in any interaction also exhibits a universal Brownian motion [13]. The main problem 
of this kind of derivation is a convincing physical origin for that universal Brownian motion, 
although a possibility is the quantum nature of space-time [14]. The chaotic behaviour of 
the solar system during its formation and evolution [15, 16] suggests a diffusion process to be 
described in terms of a Schrodinger-type equation. The description of the planetary system 
using a Schrodinger-type diffusion equation has been realized in [17]. There the authors have 
adapted the Schrodinger equation to the planetary system and shown that there exist very 
many orbits, on which possible planets may originate. That paper has stimulated us to the 
following considerations. 




QUANTIZATION AND DISCRETIZATION AT LARGE SCALES 



29 



2 Discrete distances in the gravitational field of an as- 
tronomical body 



Let us consider a body of the mass M p , which orbits a central body of the mass M and has 
the potential energy V (. x , y, z ) in its gravitational field. Because planets and moons of the 
giant planets revolve approximately in the same plane, we consider z — 0. Because they 
revolve in the same direction, we choose directions of the axes x, y and z such that the 
planets or moons of giant planets revolve counter-clockwise. Then we write the modified 
Schrodinger equation for the wave function ^ = iji(x,y) from the part of the Hilbert space 
L 2 (R 2 ) Pi C 2 (R 2 ) and the eigenvalue 0 > E G R in the form 

h 2 l d 2 d 2 \ 

-2 %{w + >wr+ nt ' v) * =E *' (5) 

where % M fa 1.48 x 10 15 M p , V(x,y) = V(x,y,z ) and E is the total energy. Negative E 
classically correspond to the elliptic Kepler orbits and the localization property (bound 
state) is conserved also in the quantum mechanics for such total energies E. The factor 
1.48 x 10 15 is not a dimensionless number, but the unit of its measurement is m 2 s -1 . With 
respect to the unusual unit we do not wonder that Agnese and Festa [5] consider this factor 
in the form of a product, such that Hm = A mcM p , where A m ~ 4.94 x 10 6 m. 

We transform equation (5) into the polar coordinates, 



till f d 2 j> 1 tty 1 d 2 ip\ 

2M p y dr 2 r dr r 2 DO 2 J 



+ V (■ r)ip = EtP, 



( 6 ) 



where iji = ip{r, 6)=i(i(r cos 9, r sind) and V (r) = l/(rcosd, rsin0) does not depend on 9. 
Particularly we choose 



V(r) = 



GM p M 

r 



(7) 



With respect to the Fourier method we assume a solution of the equation (7) in the form 



-0(r, 6) = R(r)Q(9). 



( 8 ) 



The original eigenvalue problem is transformed, equivalently, to two eigenvalue problems 



e"(0) = — a©, 

0(0) = 0(2vr) 



and 



R!'(r ) + - R'(r ) + i-4 + 



E-V{r) 2 -^ 

n M 



R{r ) = 0, 



lim [y/rR(r)\ = 0, y/rR(r) G L 2 ((0, oo)). 

i — >o+ 



The solution of the problem (9)-(10) has the form 



(9) 

( 10 ) 

( 11 ) 

( 12 ) 



9j( 0) = -s=exp (iW) 
V 2 7T 



(13) 



for l = ±V^A G Z. 

Here l = 0 should mean a body, which does not revolve at all. In the classical mechanics 
such a body moves close to a line segment ending at the central body, and it spends a short 
time in the vicinity of this body. In this paper we utilize some - not all - of the concepts 
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of quantum mechanics and we will not avoid the case l = 0 [17]. In (13) l 
corresponds to the counter-clockwise revolution. 

Respecting (7), the equation (11) becomes 

... . 1 .. ( l 2 2M„ / GM V M\ ) . , 

R (r) + -R (r) + ( — ) ) R(r) = 0, 



= 1,2, ...,oo 



(14) 



where 



Let us note that 



j-j 2M P E 
— h 2 

n M 



2 E 
(Amc ) 2 Mp 



M V GM V M GM 



ft 2 

n M 



(A mc)‘ 



On substituting r = and introducing 






equation (14) becomes 



-R W (p) + ~${p) + T 4 2 I 

P V 4 P P' ) 



where 



k = 



GM 



(A MC) 2 ^' 

For later reference let us note that, inversely, 



Expressing R(p) in the form 

we obtain an equation for u(p) 

n"(p) + 



(15) 

(16) 



(17) 



(18) 



(19) 



sfB GM 

(A Mcfk' 


(20) 


-E . (A M cf E 

Mp 2 


(21) 


(GM) 2 
2(A M c) 2 k 2 


(22) 


R(p) = u (p ) , 

Vp 


(23) 


+ --h 2 -b A1 «(p)=o, 


(24) 



where l' = l. It is familiar that this equation has two linear independent solutions Mfc //(p), 
Mk-i'(p), if l' is not an integer number. When V is integer, the solution Mk-i'(p) must be 
replaced with a more complicated solution. It can be proven that the other solution is not 
regular for p = 0 (it diverges as lnp for p — > 0). The remaining solution A4k,i(p) can be 
transformed to a wave function from the space L 2 ((0, 00 )) if and only if k — l — | = n r is any 
nonnegative integer number. We choose this function to be 



Uki(p) — CkiM kt i{p), 



(25) 
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where C k i is an appropriate normalization constant and M k j(p) is a Whittaker function, 
namely 



M k j(p) 


= p l+ ^ exp (-0 $ k + \i 21 + i;p) , 


(26) 


where <f> is the confluent (or 


degenerate) hypergeometric function. In (25) the constant C k i 


has the property 








[ r[R k i(r)fdr = 1, 
to 


(27) 


or it is 


r _ o /p 1 (n + l — 1)! 

Ckl ~ V (2l)\\2k(n-l-l)\- 


(28) 


Then 






R k i(r) = 2 VB, 


, 2fcr(n + Z) eXP ^ rvA ®)( 2rv/ ^) /L n-i- i(2rVB), 


(29) 



where n — k + L^_ l _ 1 (x) is a Laguerre polynomial, and the relation (20) holds. 



3 Interpretation of formulae derived 

Having solved the modified Schrodinger equation, we address interpretation of the formulae 
derived. The probability density Pki(r) of the revolving body occurring at the distance r 
from the central body is 

Pki(r ) = r[R k i(r )] 2 , r G [0,oo). (30) 

Mean distances of the planets are given by the relation 

/»oo 

f ki = / rP kl (r)dr (31) 

to 

= l [(2 k - n r )(2k - n r + 1) + 4n r (2A; - n r ) + n r (n r - 1)] , (32) 

where n r — n — l — 1 , k = |, |, |, . . . , oo and l — 0, 1 , 2, . . . , n. 

For the solar system M = M Sun holds and the Bohr radius of the solar system n 0 = 
0.055 AU. For survey one finds some expectation values r k i for selected values /c, l with the 
specification of described bodies in table 1 (cf. [IT]). 

Even though also in this case an empirical formula is tested for distribution of planetary 
distances, the predicted orbits fit those of the bodies in this solar system. 

Using the graphs of the probability densities we have plotted for every predicted orbit of 
this system, we obtain surprising results. The graphs of the probability densities for each 
orbit with k < | and with y < k < y are contained, respectively, in figure 1 and in figure 2. 
The vertical axis denotes the probability density P k i(r ) and the longitudinal axis designates 
the planetary distance r from the Sun. In figure 1 graph no. p = 1 is interpreted such that 
the highest probability density is assigned to the orbit of the radius of 0.055 AU and from 
the calm shape of the graph we infer that an ideal circular orbit is tested. In figure 1 graph 
no. p = 14 is interpreted such that the highest probability density is assigned to the orbit 
of the radius of 3.32 AU and, of many peaks, which wave the shape, we infer that no stable 
circular orbit is tested. After performing the analysis for all the orbits, we obtain only a 
small number of stable circular orbits. The orbits, on which big bodies - planets - may 
originate, are listed in table 2. 

It emerges that, for every number k, there exists only one stable orbit, on which a big 
body - a planet - may originate. Then we can interpret the number k as the principal 
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quantum number and l as the orbital quantum number equal to the number of possible 
orbits, but only for the greatest l there exists a stable orbit of a future body. A planet 
which does not confirm this theory is the Earth. Since the description based on the modified 
Schrodinger equation for the planetary system is not fundamental, it could not fit all the 
stable orbits. Other deviations are likely to be incurred by collisions of the bodies in early 
stages of the origin of the planets, thus nowadays we already observe elliptical orbits, which 
are very close to circular orbits. 

This procedure has been applied to moons of giant planets by us. It emerges that the 
moons of giant planets also are fitted by the modified Schrodinger equation and appropriate 
expectation values. Especially, the predicted stable circular orbits of Jupiter’s moons are 
presented in table 3. For Jupiter it holds that M = M Jup and the Bohr radius (4) of this 
system rq = 6287 km. It emerges that the predicted lunar orbits fit the measured orbits of 
the moons orbiting Jupiter. 

4 Conclusions 

In this paper we assume that there exists a law by which big objects - planets and moons of 
giant planets - do not originate anywhere, but at allowed distances from the central body. 
Unnegligible number of authors have issued from similar assumptions and derived empirical 
formulae for parameters of allowed orbits. 

The results we have presented in this paper are based on a modified Schrodinger equation, 
which has been applied to the planetary system by us for the quantum theory contained in 
the Schrodinger equation to create an interesting view of the birth of such a stellar system, 
namely the orbits of planets and moons being approximately quantized. 
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Table 1. Predicted distances of bodies from the Sun 



Body 


k 


l 


r k i [AU] 


— 


1 

2 


0 


0.055 


Mercury 


3 

2 


1 


0.332 


Mercury 


3 

2 


0 


0.387 


Venus 


5 

2 


2 


0.829 


Earth 


5 

2 


1 


0.995 


Earth 


5 

2 


0 


1.050 


Mars 


7 

2 


3 


1.548 


Hungaria 


7 

2 


2 


1.824 


Hungaria 


7 

2 


1 


1.990 


Hungaria 


7 

2 


0 


2.046 


Vesta 


9 

2 


4 


2.488 


Ceres 


9 

2 


3 


2.875 


Hygeia 


9 

2 


2 


3.151 


Camilla 


9 

2 


1 


3.317 


Camilla 


9 

2 


0 


3.372 


Jupiter 


11 

2 


0 


5.031 


— 


13 

2 


0 


7.021 


Saturn 


15 

17 

2 


0 


9.343 


Chiron 


0 


11.997 


Chiron 


19 

2 


0 


14.982 


Uranus 


21 

2 


0 


18.300 


— 


23 

2 


0 


21.948 


HA2 (1992), DW2 (1995) 


25 

2 ^ 

2 


0 


25.929 


Neptune 


0 


30.241 


— 


29 

2 


0 


34.885 


Pluto 


31 

2 


0 


39.861 
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Table 2. Bodies with stable circular orbits. 



Body 


k 


l 


ru [AU] 


— 


1 

2 


0 


0.055 


Mercury 


3 

2 


1 


0.332 


Venus 


5 

1 

2 


2 


0.83 


Mars 


3 


1.54 


Vesta 


9 

2 


4 


2.49 


Fayet comet 


11 

2 


5 


3.64 


Jupiter 


13 

2 


6 


5.03 


Neujmin comet 


15 

17 

2 


7 


6.636 


— 


8 


8.46 


Saturn 


19 

2 


9 


10.5 


— 


21 

2 


10 


12.77 


Westphal comet 


23 

2 


11 


15.26 


Pons-Brooks comet 


25 

27 

2 


12 


17.97 


Uranus 


13 


20.9 


— 


29 

2 


14 


24.055 


— 


31 

2 


15 


27.43 


Neptune 


33 

2 


16 


31.02 


— 


35 

3 ^ 

2 


17 


34.84 


Pluto 


18 


38.88 


— 


39 

2 


19 


43.134 
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Table 3. Moons of Jupiter with stable circular orbits. 



Body 


k 


/ 


r u [km 




— 


2 


0 


6287 




— 


3 

2 


1 


37722 




Halo ring 


5 

2 


2 


94305 




Outer ring 


7 

2 


3 


176036 




— 


9 

2 


4 


282915 




Io 


11 

2 


5 


414942 




Europa 


13 

2 


6 


572117 




— 


15 

17 

2 


7 


754440 




— 


8 


961911 




Ganymede 


19 

2 


9 


1.19x10 


6 


— 


21 

2 


10 


1.452x10® 


Callisto 


23 

2 


11 


1.735x10® 
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Figure 1: Probability densities for a particle in states with quantum numbers k, /, which 
correspond, respectively, (p is an ordinary number) to Mercury (p — 1, l — 1), Mercury 
(p — 2, l — 0, the second possibility), Venus (p — 3, l — 2), Earth (p — 4, l — 1), Earth 
(p — 5, l — 0, the second possibility), Mars (p — 6, l — 3), asteroid Hungaria (p — 7, l — 2), 
asteroid Hungaria (p — 8, l — 1, the second possibility), asteroid Hungaria (p = 9, l = 0, 
the third possibility), asteroid Vesta (p = 10, Z = 4), asteroid Ceres (p = 11, l = 3), asteroid 
Hygeia (p=12, / = 2), asteroid Camilla (p = 13, / = 1), and asteroid Camilla (p = 14, l = 0, 
the second possibility). Here k G . . . , the quantum number k repeats n(— k + |) 

times and r is measured in AU. 
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Figure 2: Probability densities for a particle in states with quantum numbers k, l, which 
correspond, respectively, (p is an ordinary number) to Jupiter (p — 1), nothing ( p = 2), 
Saturn (p — 3), Chiron ( p = 4), Chiron ( p = 5, the second possibility), Uranus ( p = 6), 
nothing (p — 7), HA2 (1992), DW2 (1995) (p — 8), Neptune (p — 9), nothing ( p = 10) and 
Pluto ( p = 11), k G { y, y, . . . , y), l — 0- Here k — p + |, r is measured in AU. 
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Abstract 

We propose new cosmological model of the universe on the basis of a recombination 
and an ionization of the hydrogen atom. According to the Bohr atomic model, we can 
explain main epochs in the universe in agreement with observations. We can predict 
current values of the Hubble parameter according to WMAP observation. We derive 
that the Hubble parameter can be quantized in time according to the quantum numbers 
and we express the dependence of the Hubble parameter on time. 



1 Introduction 

The connection of quantum physics with relativity is main physical problem of 21st century. 
Many authors have searched for this connection using the rescaled Planck constant [1, 2], 
Carneiro [3] has shown that there exists a rescaling factor A for the large-scale quantization, 




where T, R, and M are age, size, and mass of the universe, respectively, t, r, and m are typical 
values of time, length, and mass appearing in particle physics, respectively It emerges that 
this scaling factor is equal just to a power of ten. A scaling factor A ~ xo 38 ~ 41 is assigned to 
the size of galaxies and thus the rescaled Planck constant is H 10 81 , where H = h\ 3 . This 
is really in good agreement with the Dirac large number hypothesis. Dirac [4] has arrived at 
the revolutionary hypothesis that the gravitational constant G varies in time. Many authors 
have applied these ideas to a new cosmological model of space with better or worse success. 
Nottale [5] has discussed a cosmological constant A as sum of a general relativistic term and 
of a quantum term, scale- varying, gravitational self-energy of virtual pairs in agreement with 
present observational limits. Shemi-zadeh [6] has found the equality of fundamental large 
numbers to the exponent of the inverse value of the fine structure constant. The scaling 
law for the cosmological constant has been applied [7] in the connection with the standard 
cosmological model. Berman and Trevisan [8, 9] have derived a possible time dependence 
of some other parameters of the universe (the number of nucleons, the speed of light c, 
the gravitational constant G, and the energy density. The variation of basic fundamental 
constants has big impact on the planetary science, the cosmology, and the evolution of stars. 
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We can find very interesting description and also discussions about the time dependence of 
basic fundamental constants in [10]. 

Tifft [11, 12] has published a possible redshift quantization. He has indicated the quan- 
tization of the distances to galaxies from the Earth. Arp [13, 14] has studied the redshift 
quantization from his observations and Guthrie and Napier [15] have found the redshift peri- 
odicity in local supercluster. These results showed that the quantization at large scales could 
be possible. But other authors have critized the redshift quantization. According to [16], 
intermediate periodicities in the redshifts are due to other geometric irregularities. Hawkins 
et al. [17] have found no evidence for the redshift quantization in the 2dF survey. Bajan et 
al. [18] have published a historical review of the redshift quantization. Their conclusion has 
been that the redshift quantization is an effect which can really exist and they have studied 
a periodicity in the Local Group of galaxies and the Hercules Supercluster. 

It is clear that the question of the redshift quantization is really open and only precise 
observations could help to solve this open problem. In this paper, we would like to explain 
why the redshift quantization could exist. We propose a new cosmological model of the 
universe on the basis of the hydrogen atom. 



2 Hydrogen atomic cosmological model of the universe 

According to the standard cosmological model, we can divide a history of the universe into 
four epochs: 

• the inflation phase of the universe, 

• the time t( rec ) when the radiation separates from the matter (the recombination of 
atoms) , 

• the time t ( ; on ) of the reionization of atoms and the creation of stars of first generation, 

• the current time to with today’s value of the Hubble parameter. 

We know exact values of the age of the universe for each epoch only from some observa- 
tions. Each new cosmological model should fit these values of the observations. One main 
problem of the standard cosmological model is to determine exact value of the Hubble param- 
eter H 0 , because a correct theoretical model is missing. Without such a theoretical model, 
we are not able to predict even values of the Hubble parameter in previous stages of the 
universe. We know the value of the Hubble parameter only from the observation [19] . Based 
upon measurements of gravitational lcnsing by using the HST, a value of H 0 = 72.6 ±3.1 
km/s/Mpc has been obtained. WMAP seven-year results, also from 2010, provide an esti- 
mate of Ho = 71.0 ± 2.5 km/s/Mpc. 

We propose a new theoretical cosmological model using the properties of the hydrogen 
atom. The idea is that a history of the universe is connected with the physics at atomic scales 
(the creation of atoms, the ionization of atoms, etc.). We understand that the hydrogen atom 
is an isolated system. It means that during the evolution of the universe, the mass of the 
electron m e does not vary in time and also a speed of light c does not vary in time, because 
the small variability of these constants renders the hydrogen atomic model unstable. 

We study a simplified model of the hydrogen atom and address to the following consid- 
eration, which comes out of the Kepler second law: Areas which are swept out by the radius 
vector of planet in equal time intervals are equal, so the elementary area swept out by the 
radius vector of planet in the aphelium in the time dt is the same as the elementary area 
swept out by the radius vector of planet in the perihelium in the time dt. For the area which 
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is swept out by the radius vector of electron in the circular model, the Kepler second law is 
valid as well. It holds that 

r e v e = 2 W, (2) 

where W is the area! velocity of the electron on the Bohr orbit, r e is the Bohr radius, v e is 
an orbital speed of the electron on the Bohr radius. The relation (2) can be rewritten in the 
form 

= W (3) 



rri P rn e 



Now we will concentrate on the left-hand side of equation (3). It holds that 

2 m e 2 m e ’ 

where we substituted v e = ca e according to [20], The parameter a e = 1/137.036 is the 
fine structure constant. According to the Hubble observation, it is valid that the speed of 
galaxies depends on the distances of galaxies from an observer. The Hubble law is of the 
form 

rH 0 = v K , (5) 

where r is the distance to galaxy from the Earth, Hq is the Hubble parameter, uh is the 
Hubble speed. Multiplying equation (4) by the Hubble parameter, we get 

r e v e ca e r e , . 

-AT# = ~^^H 0 (6 

2m P 2m P V ; 



or equivalently 



r e v e H 0 _ ca e vu e 

m P 2 m P 



where Un e = r e ^r- Applying the substitution 



^e^He 



we get the parameter A in the units [ms x kg ] and it holds that 



r e v e H 0 
m P 2 



= cA 



or equivalently 



r e v e H 0 
m P A 2 



Introducing the notation r u for a distance with vh = c, 



we get the Hnbble law in the form 



2 m P A 



r u Hn = c. 



What kind of physical mechanism produces an expansion according to the Hubble law? 
We discuss equation (8). We can substitute 



Pa = 



(13) 
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where a parameter p A is in unusual units [ms -1 kg -1 ]. This parameter resembles a momen- 
tum, but with an inverse mass 1 /m e . According to this, we apply the substitution and 
transform the units in the form 

—K = m A , (14) 

m e 

where m A is the inverse mass to m e in the unit [kg] and a parameter K = 1 in the unit [kg 2 ] . 
It is valid that 

m e m A = K, (15) 

and we call K the ” inverse unit parameter” . 

We understand the inverse mass m A as the mass of dark energy, which has the repulsive 
character and now we understand the parameter P A as the momentum of dark energy in 
correct units. We can modify equation (13) in the form 



PA 



m A r C H 0 
2 K 



(16) 



where after the modification, it is valid 



Pa K 



p H ° 

Pa = m A r e — , 



(V) 



where P A is the momentum of dark energy and for the Hubble parameter Hq in the hydrogen 
atomic model it is valid that 



H 0 = 



2P A m e 
r e K ' 



(18) 



Multiplying equation (16) by a half of the Hubble parameter, Hq/ 2, we get a repulsive force 
F a of dark energy in the form 



Fa =p A K I j- 



m A r e H q 2 
4 



(19) 



3 Epoch of universe - Radiation separates from matter 



Now we can define a time f( rec ) of the universe, when atoms recombinated together. In 
another words, we understand a recombination of atoms as the recombination to stable 
atoms. The hydrogen atom is stable on the first energy level n — 1, which is equal to the 
Bohr radius of electron r e and it is valid that F e = F 01 where F 0 is the centrifugal force. 
It means that the electrical force F e between the electron and the proton is equal to the 
centrifugal force F 0 . According to this, the electron will be on the stable orbit r e . We can 
write 



F = F 

± e o 



1 e l e 2 vl 

~a T~ = me ~ ’ 

47re 0 r/ r e 



( 20 ) 



where r e is a distance between the electron and the proton, e\ is the charge of the electron, 
e\ = — e and e 2 is the charge of the proton, e 2 = e, with an elementary charge e = 1.602 x 
10~ 19 C. As the atoms recombinate when F 0 = F A , we can write equation (20) in the form 



Fe 



Fa = 




2 

(rec) 



and express the Hubble parameter in the time of recombination in the form 



( 21 ) 



H( rec) 




m A r e 



( 22 ) 
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and for the recombination time t( rec ) it is valid that 

1 

t(rec) ~ 77 

(rec) 



(23) 



where we used F c = 8.24 x 10 8 N. 
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Figure 1: Dependence of the Hubble parameter H 0 on the age T of the universe. 

The dependence of the Hubble parameter H 0 on the age T of the universe is depicted 
in Figure 1. It is obvious that the Hubble parameter decreases with the increasing age of 
the universe. Using this dependence, we can calculate main epochs in the history of the 
universe. According to these results, our model is in good agreement with the observation of 
microwave background in the framework of the standard cosmological model of the universe 
with t( rec ) = 420000 years after the Big Bang. 



4 Epoch of universe - Formation of stars of first gen- 
eration 



The ionization of atoms and the formation of stars of first generation is next main epoch of 
the universe. According to the model for the ionization of the hydrogen atom, it holds that 
n » 1 and for a distance of the electron on the orbit n to the nucleus, it holds that 
For a distance r n , it is valid that 

r n = r e n 2 (24) 

according to the Bohr model of the hydrogen atom and for F en it holds that 



F Pn = F n = Fa = 



eie 2 



47T60 r 



= nip 



m A r n HS 



(25) 



where eo is the permittivity of free space. The Hubble parameter ion ) for the ionization of 
the hydrogen atom is of the form 



H 



(ion) 







m A r n 



(26) 
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where F en = 5.15 x 10 13 N, n = 20. We get a time t(i on ) in the form 

1 

t(ion) ~ 77 • 

-tz(ion) 



(27) 
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Figure 2: Dependence of the momentum of dark energy 2 Pa on the Hnbblc parameter Hq. 

The dependence of the momentum of dark energy ‘IP a on the Hnbble parameter H 0 is 
illustrated in Figure 2. The momentum of dark energy 2Pa is linearly proportional to the 
Hubble parameter. 



0 50 l(M) 150 200 




Figure 3: Dependence of the diameter of universe r u on the integer numbers n. 

The dependence of the diameter of universe r u on the integer numbers n is depicted in 
Figure 3. The diameter of universe increases rapidly from n — 1 to n — 20 and for n > 20 
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increases slowly. For n — 1, the recombination of atoms sets in, for n = 20, the ionization 
of atoms and the formation of stars of first generation set in, for n = 178, the present times 
set in. 



0 20 40 60 KO 100 120 140 160 180 200 




Figure 4: Dependence of the Hubble parameter Hq on the integer numbers n. 

The dependence of the Hubble parameter Ho on the integer numbers n is illustrated in 
Figure 4. For n = 1, the recombination of atoms sets in, for n = 20, the ionization of atoms 
and the formation of stars of first generation set in, for n = 178, the present times set in. 
These results are in good agreement with the standard cosmological model of the universe 
with f(i on ) = 168 x 10 6 years after the Big Bang. 

5 Current value of the Hubble parameter 

For a prediction of current value of the Hubble parameter, it is neccessary to find a correct 
mathematical relation. For an illustration, we can use n » 1 in equation (26). According to 
the observational WMAP data, the value of the Hubble parameter should be H 0 = 71.0 ±2.5 
km/s/Mpc. Applying this tolerance to the graph, we can see a connection of the momentum 
of dark energy 2Pa with the Hubble parameter Hq. Five possible results fall into the area 
of tolerance for the current Hubble parameter according to WMAP data (see Figure 5 and 
Table 1). 



The prediction of current value of the Hubble parameter H 0 is plotted in Figure 5. 
On the graph, we can follow the tolerance for the Hubble parameter according to WMAP 
observation. Five possible results fall into this tolerance for n = 178, 179, 180, 181, 182. 
After the modification of equation (16), we get 



PA 



m A r n Ho 
2 K 



and for the Hubble parameter Hq it holds that 



_ 2Kp A 
m A r n ' 



(28) 



( 29 ) 
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Figure 5: Prediction of current value of the Hubble parameter Hq. 



Table 1: Predictions of the Hubble parameter 



n 


H 0 [s' 1 ] 


Pa 


2P a 


Hq [km/s/Mpc] 


178 


2.37 x 10~ 18 


69.05 


138.11 


73.36 


179 


2.35 x 10~ 18 


68.28 


136.57 


72.54 


180 


2.32 x 10~ 18 


67.53 


135.06 


71.74 


181 


2.29 x 10“ 18 


66.78 


133.57 


70.95 


182 


2.27 x 10“ 18 


66.05 


132.11 


70.17 
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The comparison of results obtained in the framework of the standard model of the universe 
with those obtained in the framework of the hydrogen atomic cosmological model of the 
universe can be performed on the basis of Tables 2 and 3. 



Table 2: Standard cosmological model 



Parameter 


Marking 


Value 


Comments 


Hubble’s parameter 


FI o [km/s/Mpc] 


71±2.5 


only from measurements 


Age of universe 


T u [years] 


13.4T0.3 x 10 9 


calculation 


Recombination of atoms 


t (rec) [years] 




380000 


from observation 


Ionization of atoms 


t ( ion) [years] 




200 x 10 6 


from observation 


Total mass 


Aftot 


1±0.02 


including dark matter 



Table 3: Hydrogen atomic cosmological model 



Parameter 


Marking 


Value 


Comments 


Hubble’s parameter 


H 0 [km/s/Mpc] 


72.54 


calculation for n = 179 


Age of universe 


T u [years] 


13.44 xlO 9 


calculation 


Recombination of atoms 


t( rec) [years] 




420978 


calculation for n — 1 


Ionization of atoms 


t(i on) [years] 




168 x 10 6 


calculation for n = 20 


Total mass 


Af-tot 


1 


product of dark energy and mass 



6 Epoch of universe - phase of inflation 

In our model, we define a phase of inflation if F A » F e . In this phase, the repulsive 
force Fa is much higher than the Coulomb force F e and therefore atoms could not be stable 
and exist. It is very difficult to predict the value of the Hubble parameter in the epoch of 
inflation. It is not possible to apply the hydrogen atomic cosmological model. However, 
we understand that the calculation of the Hubble parameter in this epoch can be based on 
a model of nucleus of the hydrogen atom with nuclear forces. According to the hydrogen 
atomic cosmological model, the speed of light does not vary in time during all epochs, because 
of the recombination of atoms. 



7 Conclusion 

We propose new theoretical cosmological model of the universe in the analogy with the 
hydrogen atom. According to the Bohr model of the hydrogen atom, we characterize main 
epochs in the universe for the time of the recombination of atoms (t( re c)) and that of the 
reionization of atoms f( ion )- We predict also current values of the Hubble parameter according 
to WMAP measurements. We derive the relation for temporal dependence of the Hubble 
parameter. In this model, we discuss also properties of dark energy. According to our study, 
the Hubble parameter does not change continuously in time, but it changes discontinuously 
in the dependence on the quantum numbers n. We provide the comparison with results 
obtained in the framework of the standard model of the universe. 
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Abstract 

The recent experimental findings have shown that our understanding of the solar system is 
surprisingly fragmentary. As a matter fact, so fragmentary that even new physics might find place 
in the description of phenomena like the precession of equinoxes and the recent discoveries about 
the bullet like shape of heliosphere and strong magnetic fields near its boundary bringing in mind 
incompressible fluid flow around obstacle. TGD inspired model is based on the heuristic idea that 
stars are like pearls in a necklace defined by long magnetic flux tubes carrying dark matter and 
strong magnetic field responsible for dark energy and possibly accompanied by the analog of solar 
wind. Heliosphere would be like bubble in the flow defined by the magnetic field inside the flux 
tube inducing its local thickening. A possible interpretation is as a bubble of ordinary and dark 
matter in the flux tube containing dark energy. This would provide a beautiful overall view about 
the emergence of stars and their helio-spheres as a phase transition transforming dark energy to 
dark and visible matter. Among other things the magnetic walls surrounding the solar system 
would shield the solar system from cosmic rays. The model leads to a vision about formations of 
stars and galixes as ’’boiling” of dark energy to matter. Also a model for the cosmic rays emerges 
allowing to identify the acceleration mechanism using recent findings about cosmic rays. 
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1 Motivations 

The inspiration to this little contribution came from a discussion with my friend Pertti Karkkainen 
who told me about the work of Walter Cruttenden m- Cruttenden is a free researcher working with 
an old problem related to the astronomy of the solar system, namely the precession of equinoxes [9], 
Equinoxes [2J correspond to the two points at the orbit of Earth at which the Sun is in the plane 
of the equator (if Earth’s spin axes were not tilted this would be the case always). What has been 
observed is an apparent movement of fixed stars relative to the Earth bound observer. The period of 
the equinox precession is about 26,000 years. The angular radius of the precession cone is about 23.5 
degrees. The rate of precession is approximately 50 arc seconds per year but is not strictly constant. 

The precession of equinoxes reduces to precession which is a well-known phenomenon associated 
with the motion of a rigid body with one point fixed. Precession [8] means that the spin axis of the 
spinning system rotates around fixed axis along the surface of a cone. One can distinguish between 
a torque free precession and precession induced by torque. Precession can be accompanied by a 
nutation [6] : the tilt angle of the spin axes with respect to fixed axes varies with time. The nutation 
for Earth is well-understood process determined by the local gravitational physics. In the case of 
precession the situation is not so clear. 

1.1 Two basic theories explaining the precession of equinoxes 

There are two basic theories of precession. 

1. The precession of equinoxes could be governed by a local dynamics being due to the precession 
of the Earth with respect to solar system. Earth is indeed a prolate ellipsoid and the precession 
would be caused mainly by the gravitational fields of Sun and Moon (lunisolar model). According 
to the summary of Cruttenden HU, Newton’s equations did not work and d’Alembert and others 
have added and changed input values to fit the observed precession. The latest 2000A version 
includes almost 1400 terms but it still fails to accurately predict variations in the precession 
rate. The theory is also plagued by a ’’measurement paradox”. Studies show that the changes in 
Earth’s orientation relative to Sun and other planets are small (few arc seconds per year instead 
of 50 arc seconds) as compared to the equinox precession. 

2. The precession of equinoxes could be also due to the precession of the entire solar system 
regarded as a rigid body with one point fixed and would be caused by some hypothetical binary 
companion of Sun. Usually the binary companion is thought to be star of planet like system but 
this is not necessary. This model is known as binary model and was first proposed by Indian 
astronomer Sri Yukteswar, The predicted period was 24,000 years. According to the summary 
of Cruttenden , the binary model of Yukteswar has turned out to be more accurate over 100 
year period Kl- 
in principle the observation of the precession from some other planet could select between the two 

approaches. If the precession were similar at two planets then the precession of the entire solar system 
would be strongly favored as an explanation of the equinox precession. 

1.2 Some hints 

The basic challenge for the binary theory is of course the identification of the binary. There are some 
hints in this respect listed by W. Cruttenden in the articles at his homepage. Consider first what has 
been learned from the structure of heliosphere during last years. 

1. The data from Voyager 1 and Voyager 2 have revealed that heliosphere is asymmetric jT6]. The 
edge of the heliosphere (the place where the solar wind slows down to sub-sonic speeds and is 
heated) appears to be 1.2 billion kilometers shorter on the south side of the solar system than 
it is on the edge of the planetary plane. This indicates the heliosphere is not a sphere but has 
a shape of a bullet. In a sharp contrast with the naive expectations, the magnetosphere of Sun 
would not be like that of Earth which is compressed on the day side by solar wind and has a 
long tail on the night side. 
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2. There is also evidence from Voyager 2 for a strong magnetic field [TS]. Also the temperature just 
outside the boundary zone defining the boundary of the solar inner magnetosphere was ten times 
cooler than expected. The presence of the strong magnetic field is not easy to understand since 
the interstellar space consists of extremely tenuous gas. The proposal is that the interstellar 
magnetic field could be forced to flow around the helio-magnetosphere much like fluid flows 
around obstacle. This increases the density of flux lines and interstellar magnetic field would 
become stronger locally. Heliosphere would be like a bubble inside magnetic flux tube expanding 
it locally. 

The direction of the local magnetic field at the edge of the heliosphere differs considerably from 
that for the interstellar magnetic field thought to be parallel to the galactic plane. The tilt angle 
is about 60 degrees. Therefore one can challenge the identification of the strong local magnetic 
field as galactic magnetic field. 

3. Between June and October 2007, the STEREO spacecraft '\V1\ ’’detected atoms originating 
from the same spot in the sky: the shock front and the helio-sheath beyond, where the sun 
plunges through the interstellar medium, and found energetic neutral particles from beyond the 
heliosphere that are moving towards the sun M- This would suggest magnetic flux tube like 
structure and the flow of neutral particles along the flux tube towards the Sun so that an analog 
of solar wind would be in question. 

Also the behavior of comets suggests that the understanding of the solar system is far from com- 
plete. The behavior of the comet Sedna thought to belong to the inner Oort cloud ,7] cannot be 
explained in terms of theory assuming only solar and planetary gravitational fields. Typically comets 
move along periodic orbits returning repeatedly near some planet of solar system (typically Neptune) 
which has kicked the comet to its highly eccentric orbit. Sedna PQ (thought to be a ’’dwarf planet”) 
seems to be an exception in this respect. Sedna has an exceptionally long and elongated orbit (aphe- 
lion about 937 AU and perihelion about 89.6 AU), period is estimated to be 11,400 years, and Sedna 
does not return near any planet periodically as the assumption that it belongs to the scattered disk 
would require. 

What could be the origin of Sedna? 

1. It has been suggested that that Sun has an dim binary companion - christened Nemesis [5] - 
at a distance of thousands of AUs. This companion could explain the behavior of Sedna, and 
has been also proposed to be responsible for the conjectured periodicity of mass extinctions, the 
lunar impact record, and the common orbital elements of a number of long period comets. 

2. Second proposal is that Sedna has been kicked to its orbit by some object. This object could be 
an unseen planet much beyond the Kuiper belt [3] (Kuiper belt is outside planet Neptune and 
extends from 30 AU to 55 AU). It would have mass about 5 times the mass of Jupiter and be at 
distance of roughly 7850 AU from the Sun in the inner Oort cloud. It could be a single passing 
star or one of the young stars embedded with the Sun in the stellar cluster in which it formed. 
This might have happened already in the Sun’s birth cluster (cluster of stars). 

3. Also the behavior of the comets in outer Oort cloud (very eccentric orbits and long orbital 
periods) might reflect the influence of a binary companion whose mass distribution is such that 
this kind of orbits are generic. For spherical objects one would expect nearly circular orbits. 
String like object would satisfy this condition as will be found. 

1.3 The identification of the companion of the Sun in the framework of 
standard physics 

Consider first the identification of the companion of the Sun responsible for the precession of the solar 
system as a whole but staying in the framework of the standard physics. In this context only objects 
with a spherical symmetry can be considered. 

1. The strange behavior of Sedna suggests that binary could be an unseen planet at distance of 
about 7850 AU in the inner Oort cloud. Note that Oort could extend up to 50,000 AUs which 
corresponds to .75 ly whereas the closest star - Proxima Centauri- is at distance of about 4.2 
light years. 
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2. The identification of the binary as the hypothetical Nemesis might explain the analog of the 
solar wind. If the dim Nemesis is at the same distances as the hypothetical planet, its mass 
would be only .5 per cent of solar mass. 

3. An analog of solar wind flowing along magnetic flux tubes could also come from some other star, 
say Proxima Centauri [10] . Proxima Centauri is however too light as red dwarf and too distant 
to induce the precession of the solar system as whole. 

1.4 The identification of the companion of the Sun in TGD framework 

In TGD framework one can consider more speculative ideas concerning the identification of the binary 
of the Sun. 

1. In TGD Universe dark matter and dark energy can be understood as phases of matter with large 
Planck constant [1]. For the dark energy assignable to the flux tubes mediating gravitational 
interaction between Sun and given planet the value of the Planck constant is of order GMm/v o, 
where v$/c ~ 2 -11 holds true for the inner planets. For dark matter the value of Planck constant 
is much smaller integer multiple of its minimal value identified as the ordinary Planck constant. 
Whether only magnetic energy should be counted as dark energy or whether also dark particles 
with a gigantic value of Planck constant should be identified as dark energy is not quite clear. 

2. Magnetic flux tubes are identified as carriers of dark matter. This hypothesis plays a key role 
in TGD inspired quantum biology and cosmology. The flux tubes can have arbitrary large 
length scales. During the cosmology space-time would have consisted of cosmic strings of form 
X 2 x Y 2 C M 4 x CP2 with X 2 minimal surface and Y 2 complex sub-manifold of CP-2- In 
the course of the cosmic evolution their M 4 projection would have become 4-dimensional and 
they would have become magnetic flux tubes. The proposal is that galaxies are like pearls in a 
necklace formed by flux tubes p]. 

The density pdark of the magnetic energy is enormous for cosmic strings: the length L of cosmic 
string corresponds to a mass which is a fraction G/HqR 2 ~ 10 -4 of the mass of a black hole 
with radius L. The thickening of the cosmic string to a flux tube respects the conservation of 
the magnetic flux so that the strength of the magnetic field scales down like B oc 1/S, where S 
is the area for the transversal cross section of the flux tube. By a simple scaling argument the 
density of the magnetic energy per unit length of the flux tube scales down like dE m /dl oc 1/5. 

If energy is conserved if the length of the cosmic string scales up like 5 in the cosmic expansion: 
d oc y/L proportionality analogous to that encountered in the case of diffusion would relate to 
each other flux tube radius and length. Also the primary p-adic length scales L p assignable to 
particles and the secondary p-adic length scales L p> 2 characterizing the corresponding causal dia- 
mond CD relate in a similar manner. This would suggest that the p-adic length scale assignable 
to a given particle (of order Compton length) corresponds to the thickness of the magnetic flux 
tube(s) assignable to the particle and the size of CD to the length of the(se) magnetic flux 
tube(s). Similar scaling holds true for the density of dark matter per unit length of the flux 
tube. 

The dark matter associated with the flux tubes would generate transversal 1/p gravitational 
field explaining the constant velocity spectrum of distance stars in the galactic halo. The basic 
prediction is free motion along the direction of the cosmic string perturbed only by the mass of 
the galaxy itself. 

3. The fractality of the TGD Universe suggests the pearls in the necklace model applies also to 
stars. The magnetic flux tube idealizable as a straight string would be roughly orthogonal to 
the plane of the planetary system possibly associated with the star and the spin axis of the 
star would be nearly parallel to the flux tube. If one combines this picture with the previous 
discussion, the simplest proposal is obvious. The binary companion of the Sun is the magnetic 
flux tube containing dark matter. 

Newtonian theory for the gravitation in planetary system works excellently and this poses strong 
constraints on the pearls in a necklace model will be discussed in more detail. 



